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Abstract 

Electroweak radiative corrections can be evaluated in the Sudakov approximation, a systematic 
high energy expansion known to be relevant for the analysis of future collider experiments in 
the TeV energy range. In the Minimal Supersymmetric Standard Model and at next-to-leading 
order, Sudakov electroweak corrections satisfy remarkable relations at the one loop level. Explicit 
computations in component fields are available for various different 2 — ► 2 processes relevant for 
Linear Collider or LHC physics. The Sudakov corrections turn out to be equal or closely related 
in several classes of processes differing by the replacement of certain final or initial states with 
their superpartners. This fact suggests that supersymmetry is partially restored at high-energy. 
We analyze the supersymmetric structure of such relations by computing the Sudakov corrections 
in the framework of superfield perturbation theory. As a simple application, we derive in full 
details an extended complete set of supersymmetric relations among different processes related by 
supersymmetry to the fundamental fermion pair production process e + e~ — > ff. 

PACS numbers: 11.30.Pb, 12.15.Lk, 11.15.Bt 
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I. INTRODUCTION 



The calculation of radiative corrections in models of fundamental interactions will play 
an important role in the theoretical analysis of future collider experiments. In this paper, 
we work in the Minimal Supersymmetric Standard Model (MSSM), where higher order 
perturbative effects can be calculated systematically due to renormalizability. 

A general discussion of the structure of radiative corrections in the MSSM can be found 
in Q. Actua. captions involve the e X cha„ g e of virtua, state, including of course, the 
various supersymmetric partners of Standard Model particles. In principle, supersymmetry 
(SUSY) could simplify the analysis by providing an underlying structure. In practice, su- 
persymmetry is softly broken at low energy making very hard to exploit the symmetry in 
the actual evaluation of definite physical processes |2|. 

Thus, most of the existing calculations of radiative corrections in the MSSM are performed 
in the familiar component formalism. Physical fields are treated separately and the SUSY 
multiplet structure is not fully exploited in the perturbative expansion. 

This is an unpleasant situation, but the following considerations suggest the possibility 
of recovering supersymmetry to some extent in the phenomenological analysis of radiative 
corrections. In the future collider experiments in the TeV energy range, in particular the 
Large Hadron Collider (LHC) and International Linear Collider (ILC), the typical experi- 
mental situations will be characterized by large invariant masses. A high energy expansion 
seems therefore a natural opportunity to achieve major simplifications in the structure of 
loop corrections. 

Indeed, high energy expansions enjoy several useful features in this respect. First, helicity 
conservation reduce the number of relevant amplitudes. Second, it is known that radiative 
effects beyond production thresholds are rather smooth and can be described by the so- 
called logarithmic Sudakov expansion (LSE) [3]. In particular, in a MSSM scenario with 



light SUSY partners below the typical value of 300-400 GeV, such expansions effectively 
apply at energies in the TeV range. The practical consequences of LSE have been deeply 
analyzed in the Standard Model § and in the MSSM both at ILC Q, H, E | and at LHC fl. 
Radiative corrections are expanded in a logarithmic variable that, in the electroweak case, is 
typically log(s/M^), where s is the squared center of mass energy and aw = ot/ sin 2 9w 
with 9w being the Weinberg weak mixing angle. Logarithms arise from ultraviolet and mass 
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singularities. At one loop and next-to-leading order (NLO), there appear contributions of 
the form aw log 2 (s/M^) and aw^og{s/ M w ) which have been shown to be sizable at the 
planned energies. For instance, at y^i = 1 TeV they give typical corrections of tens per cent 
to the one loop cross sections. 

The accuracy of LSE has been discussed by explicit comparisons with full one loop calcula- 
tions of specific processes In practice, LSE leads to very simple and compact expressions 
that allow rather straightforward phenomenological analysis [9]. With a conservative atti- 
tude, processes where LSE predicts large observable effects have strong motivations for a 
full treatment. 

We now come back to the issue of supersymmetry and consider what happens in the 
MSSM. Indeed, in this case we observe quite remarkable special features of LSE. The most 
important is that LSE turns out to be independent from most of the soft breaking terms at 
NLO. Therefore it isolates the dependence of observables on very specific parameters, like 
for instance the vacuum alignment angle tan (3 j^| . 

Independence from soft breaking can be argued by dimensional arguments in terms of the 
high energy mass suppression of low dimensional breaking operators. It suggests that some 
amount of SUSY could be recovered at the level of LSE. This conjecture has been checked 
as a by-product of complete one loop calculations of various processes relevant for ILC and 
LHC |7[ in the component approach. 

The interesting result is that Sudakov corrections to processes involving particles in the 
same SUSY multiplets are closely related, as explained in full details below. We shall denote 
briefly such kind of results Sudakov Supersymmetric Relations (SSR). In this paper, we shall 
discuss a proof of such relations within the framework of superfield perturbation theory. 

The plan of the paper is the following. In Sec. |H] we discuss the specific details of 
SSR. It will emerge that the Sudakov correction can be split into three contributions, so 
called universal, Yukawa and angular dependent. In Sec. IIHI we illustrate the supergraph 
calculation of the universal correction. In Sec . 11 VI we discuss the Yukawa correction. Finally, 
in Sec. Elwe illustrate the evaluation of box supergraphs that compute the angular dependent 
correction. In Sec. IVIl we summarize our results and discuss future developments. 
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II. SUDAKOV SUPERSYMMETRIC RELATIONS 



We now review the observed SSR as they emerge from explicit one loop calculations in 

Q 

the MSSM in component fields. We closely follow the notation of |7( and work at one loop. 
We begin the discussion by considering the ILC process (/ 7^ e) 

e«e" -> fpjp, (1) 

where a,/3 are chirality left (L) or right (R) indices. We denote by s the squared center of 
mass (cm.) energy of the initial e + e~ pair. 

The electroweak Sudakov correction to the amplitude for the process (JH) is described at 
high energy and at next-to-leading logarithmic order in the following very simple form 

^Sudakov = A Bor n(1 + ^ + TJ + V + ^ + ^RC (2) 



where the initial (c^) and final (c¥) universal coefficient c u is 



/3 ~ /3 ' u a/3 



where J, T 3 are the total and third component of isospin and K is the hypercharge. We 
denote by g and g' the SU(2) x U(l)y gauge coupling, following standard notation. 

As we claimed, Eq. (J2J) is valid at NLO, i.e. when we neglect all terms which are not 
growing with the cm. energy. The mass scale of the linear logarithm is not fixed at NLO in 
the logarithmic expansion. Instead, the mass scale in the squared logarithm is not arbitrary. 
In components, several diagrams with exchange of gauge bosons do contribute. The mass My 
can be M 7 , Mz or My/ according to the exchanged virtual gauge boson. For the photon, we 
set fictitiously M 7 = Mz- This permits a unified treatment of the neutral current electroweak 
interactions. In the application to a physical process, the correct procedure will require to 
reinsert M 7 back into inclusive observables, include initial and final radiation effects, and 



finally take the physical M 7 — >• limit 
scale separation at two loop level in 



see for instance the detailed discussion about mass 

m 

The two terms in (jHJ) have a clear separate origin from the neutral (7, Z) and charged 
(VP A± ) current sectors, as shown by the Casimir operator of U(l)y or SU(2). When the 
approximation Mz — Mw is used, one recovers the above simple formula where the isospin 
and hypercharge contributions are collect together. If M z and M w are kept distinct, the 
two contributions are separated accordingly p. 



For simplicity, in this paper, we focus on the simpler neutral currents sector to illus- 
trate the superfield technique. More precisely, we aim at recovering the SSR involving the 
hypercharge U(l)y part of the gauge group, i.e. for c u , 

- tL. 9 2 ¥ (2 log * - log 2 °) . (4) 



16vr 2 y 4 V M l M h 
Intrinsic charged currents contributions can also be considered, but there are minor compli- 
cations due to the triple gauge boson vertex. The Sudakov correction is gauge invariant, and 
it is convenient to organize the various diagrams in a gauge invariant way by means of the 
pinch technique jllj]. We leave details regarding the implementation of the pinch technique 
with superfields to a separate work. 

The Yukawa coefficient c Y is present if / is heavy and we keep its mass fljj]. Let us focus 
on the important case of top or bottom quark production. If we denote by /' the SU(2) 
partner of /, we have 

Y _ 9 2 ( 1 + 6 P , R mj 5 PtL m}, \ s 
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where rh t = m t / sin/3 and = mj/ cos (3, and M is a mass scale not fixed at NLO. It is 
typically taken as M = m t . It is important to remark that it is only and precisely from this 
correction that a dependence on the MSSM parameters arises at very high energy in this 
process. 

The angular dependent correction is 

t 



c an S = —\ofr S 



afB ' 16tt 2 6 M 2 



(tan9^Y a Y p + AT^) log-+ 

Va,L0p,L I r , —t r , —U 
1 d b j log d t j log 



(6) 



1/4 tan 6 2 w Y a Yp + T3T* \ °' J b s J 

where t and u are standard Mandelstam variables. The hypercharge contribution is partic- 
ularly simple 

'2 y Y a f 

Crr ~ 4tt 2 2 2 l0g Mf [ ° g u {n 

Finally, the Renormalization Group (RG) contribution A RG contains the logarithms due 
to RG running of the coupling constants. This piece will be considered as known and not 
discussed in the following. It can be found in details in []|. 

These expressions are remarkably simple and can be used as a starting point for the 
analysis of the phenomenological consequences of radiative corrections, at least in the high 



energy regime. A detailed discussion of their accuracy, resummation properties, and rele- 
vance to the MSSM parameter space analysis can be found in the papers we have cited in 

fl 

the Introduction (see in particular the review j^]). 

In this work, we shall be concerned with the SUSY structure of these corrections. Such a 
structure is illustrated by the following results. Instead of (JTJ), we could consider the related 
process of sfermion production, again at ILC, 

e + a e~ - ] J;. (8) 

Here, the index (3 is not a chirality index but just denotes fp as the sfermion partner of the 
chiral fermion fp. The remarkable result is that the Sudakov correction takes the same form 
as for fermion production with the same coefficients, both for the universal, Yukawa, and 
angular dependent parts. 

Also, if we are interested in charged Higgs production e+e~ — > H + H~, we find again 
the same set of corrections under the assignment Q(H~) = —1, T 3 (H~) = —1/2. The only 
change is in the Yukawa term that reads instead 

= - 3 ■ \m cot2 v + 4; w n iog w < 9 > 

Similar results are obtained for final charginos, neutralinos 19], as well as for the initial states 
which appear in LHC processes, i.e., quarks, gluons |8j]. In the following, we shall denote by 
SSR (Sudakov SUSY relations) the equality of Sudakov corrections for processes that differ 
by the replacement of some particle in the initial or final state with its superpartner. 

The very existence of SSR proves that some amount of SUSY is recovered at the level of 
Sudakov correction, although supersymmetry is softly broken in the MSSM. As we said in 
the Introduction, we can argue that at high energy soft breaking terms are irrelevant and 
mass suppressed and SUSY relations arc recovered. An application of this idea can be found 
in [3| where the structure of gauge boson helicity conservation in the MSSM is analyzed 
under such hypothesis. 

In principle, under the assumption of unbroken SUSY, some partial information can be 
derived from Ward identities. For instance, to match the corrections to the production of a 
fermion pair to those in the sfermion case, one could start from the relation 

<W (0|e+e-^|0) = 0. (10) 
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The variation of the final state (fpfp) connects the amplitudes for the two processes e+e~ — > 
fpfp and e+e~ — ► fpfp- However, the terms coming from the variation of the initial state 
(e+e~) also contribute due to the possibility of exchanging a neutralino in the s channel. 
Thus, Ward identities suggest that SSR can be extended to include the gaugino form factor 
A — > ip*ip in the equality between the Sudakov corrections to the vector form factors 

— > ipijj and —>■ <£>*y?. This will be immediately proved in the later calculations. Beside 
these considerations, it must be checked in some details that soft terms are actually irrelevant 
and do not contribute at the Sudakov level to the right hand side of the variation Eq. (|10j). 

It is clear that that to obtain more detailed results on the Sudakov logarithmic corrections 
we should exploit in finer details their specific kinematical origin. Let us discuss the infrared 
contributions that, for instance, lead to squared logarithms. In spontaneously broken gauge 
theories the gauge boson masses My are natural infrared regulators. Thus, at large energy 
E ^> My, we can expect mass singularities to show up as powers of log(E/My), i.e. 
the infrared part of Sudakov corrections. In the Standard Model, it is well known how to 
extract these (logarithmic) singular terms by a careful analysis of the relevant diagrams Q| . 
Besides, the lowest order results can be improved by RG techniques in order to resum higher 
order effects and prove factorization properties. Advanced examples of these resummation 
techniques include the treatment of angular dependent corrections whose numerical relevance 
in actual experiments is well known to be non negligible j^, and that, as we reviewed, 
also satisfy supersymmetric relations. In principle, the extension to MSSM could be done 
in order to study the possible validity of SSR beyond one-loop. 

Here, we begin with the one-loop contributions and do not consider resummation issues. 
Instead, since we are interested in SSR, we show how the mass singularities can be isolated 
in a manifest supersymmetric way. We stick to a detailed diagrammatic calculation because 
we believe that this is the clearest framework where soft breaking independence arises in a 
controlled and explicit way. However, we do not work in components, but exploit instead 



superfield perturbation theory jlo. Il7|. 



This technique is not much used in practical calculations for phenomenological studies in 
the MSSM precisely because of the large amount of SUSY breaking. In principle, soft terms 
can be included by the spurion technique . The method is somewhat cumbersome apart 
from the analysis of UV divergences or to the (related) construction of Renormalization 
Group evolution equations [19;]. For applications, component field calculations are certainly 
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preferred. 

Instead, Sudakov corrections can be computed safely and efficiently in the superfield 
approach, as this paper will show. Indeed, the calculation is also largely simplified for the 
following reasons 

1. Scheme independence: at one loop, Sudakov corrections are physical and renormal- 
ization scheme independent. Any scheme dependence drops out because scheme de- 
pendent subtractions do not introduce contributions growing with the energy. All 
the subtleties related to regularization consistency and to the need of restoring Ward 
identities 2£| do not show up at the next-to-leading logarithmic level. 

2. Mass independence: we can insert fictitious masses in the superdiagrams, at least at 
NLO. In all diagrams leading to Sudakov corrections, the mass terms are suppressed 
at high energy as a consequences of elementary properties of the basic loop integrals. 
The only relevant mass scales entering the expressions for Sudakov corrections are the 
gauge boson masses unaffected by soft terms. 

We now consider the various supergraphs contributing the Sudakov corrections of various 
kinds. For each case, we compute the Sudakov expansion of the one-loop effective action in 
a manifest supersymmetric way and show how SSR are recovered quite easily. 

III. UNIVERSAL SUDAKOV CORRECTION 

We begin with the 1/$$ term in the one-loop effective action T. The supergraph giving 
the Sudakov correction is shown in Fig. for the final state correction. Of course, the 
same (crossed) diagram give also the correction to the initial state. We want to retrieve the 
following result 

r Su_dakov = C U ptae^ (n) 

where c u is given in Eq.®. After expansion of the superfields V, $, and this proves SSR 
for the gauge universal Sudakov correction. 

The value of the supergraph is (see for instance the analogous computation in 0(4) 
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O(q) 




3>(p) 



FIG. 1: Diagram giving the universal Sudakov correction to the V&<& form factor, 
extended super Yang-Mills |1 T] ) : 

r- - f d4k 1 ( i2) 

w " 2 / (27rY D 1 D 2 D 3 [ ] 



I 



(TO <S>(p,9) $(q, 



(k + q) 2 -\D(f + ti)D-^D 2 rf 



V(r,9), 



2 w n/ 16 

where g is the neutral coupling and D\ = k 2 — M 2 z , D 2 = (k+q) 2 — M|, _D 3 = (k— p) 2 — M|. 
The convention for chiral and vector superfields is 

$ = lf + 0^ + 0^F- iOa^e - % -B 2 Qo» - -Op 6 2 f, (13) 
V wz = 0a»6 - i6 2 6\ + itfe\ + -D 2 tf , (14) 

where, for simplicity, we have written V in the Wess-Zumino gauge. This is possible for 
external V fields in gauge invariant amplitudes. Instead, virtual V exchange will involve all 
the unphysical degrees of freedom of V. The covariant derivatives D, D acting on V are 
conveniently written after contraction with auxiliary spinors (, ( and read (if = cr M r^) 

(D = (d + CP, W = -dC - OH- (15) 

We now analyze the three terms inside the square brackets in T v ^. The first term is 
proportional to the Passarino-Veltman (PV) function B evaluated at one of the external $ 
squared momenta. It is divergent and must be regulated but does not give Sudakov loga- 
rithms being independent from the cm. energy. Basic definitions and asymptotic expansions 
of the PV functions can be found in full generality in 21]. 



In the second term of we require the finite loop integral (in the large r 2 3> M 2 ,p 2 , q 2 
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regime) 

2it) (k 2 - M 2 )[(k + q) 2 - M 2 }[(k + q + r) 2 - M|] ~ 16tt 2 1 + 12Ph ' 1 ' 
The PV functions Cn and G\ 2 have the following asymptotic expansion where ~ denotes 
equality modulo terms not growing with r 2 



1 / _ , r 



2 



Cii - 53 ( 2 log ^-V^), (17) 



2r 2 \ My My 



G> ~ -pr (18) 



The third term can be integrated by parts moving the D and D operators on the chiral and 
antichiral fields $, $. By chirality, we have D a § = = 0. Also, the expressions D 2 & and 
D 2 & are not vanishing but are mass suppressed as a consequence of the equations of motion 
of $, $. To see this, we contract the effective action vertex V$$ with the propagator 

(o|t{$( p aa) $(- P ,e 2 ,e 2 )} |o> = ^-^tm 2 d\d\5 12) (19) 

where 8 12 = (6 1 - 9 2 ) 2 (9 1 - 6 2 ) 2 . If we now act with D 2 and exploit D 2 lf D 2 = 16DD 2 , we 
see that, apart from mass corrections, any effective vertex with a D 2 Q or D $ term has not 
the necessary pole to survive the on-shell projection. 

The Sudakov correction comes therefore entirely from the second term 

rjgf^w) = z -fJ -0^ Jd'9 HqMp) (k + q ) a6l D«D & V(r). (20) 

We perform the loop integral and replace k + q — > (— i/lQii 2 )(Ciiq + C\ 2 p). Integrating by 
parts the D, D operators we find 



r^? ov (g,P,r)= (21) 

9a f J4/ 



D $(g) D a $(p) + $(q)D D a $(p) (C^ + C 12 f>)aaV(r). 



4 • 167T 2 

Since we are interested in the Sudakov enhanced terms, we can further simplify. Indeed, we 
observe that the combination (fy^d^D^ and the similar one for $ are mass suppressed as 
follows again from Eq. (|19)l . The first term of Eq. (J2T| does not contribute. Anticommuting 
D and D in the second term and taking the trace we obtain 

itS? "W. = if|a(9 • P) C'u / d ^ $ W W (22) 
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We can now replace C\\ by its asymptotic expansion and identify go = 2(g'Y/2) where g'Y/2 
is the U(1)y standard coupling. A factor 2 is needed to conform to standard normalization 
of the gauge field in components. The final result is 

r^°- = ijp (V j) 2 (2 log -J - log 2 £)/«**v-* iVU- (23) 

We have obtained in a very straightforward way the SSR for the neutral part of the correction 
c u . The scale My has been set to M z since the possible soft corrections to the gaugino mass 
cannot affect the squared logarithm scale. As a technical comment, the unphysical scalar 
component of V does propagate along the internal V line, but it can be checked that it 
does not give Sudakov logarithms. Thus, the result in Eq. (|2^j) is the correct physical gauge 
invariant one. 

To summarize, the superfield calculation of the effective action has taken into account 
simultaneously the correction to the three specific vertices appearing in the expansion of 
/ d 4 6 V wz $$, i.e. 

The first two appear in the initial and final universal correction to the amplitudes for pro- 
ducing a // or //* pair. They have been proved to be equal in agreement with SSR. The 
third one is a byproduct of the superfield calculation and, by supersymmetry, gets an iden- 
tical correction. Hence, the SSR can be extended in supersymmetric way to the gaugino 
form factor A — > <f)*ip in full consistency with the Ward identity (jl(Jj) . As a check, an explicit 
calculation in components confirms this result. 



IV. YUKAWA SUDAKOV CORRECTION 



The Yukawa corrections proportional to the squared heavy quark masses arise from a 
different kind of supergraph and have special features that we now discuss in details. We 
start from the corrections ~ m|. Later, we shall include those ~ m\. The Sudakov Yukawa 
corrections can be treated in a simplified and independent sector of the MSSM governed by 
the Lagrangian 

L = J2 I d 4 B$ e 9oV $ + y t e al3 I d 2 6 (u Q a H u p + h.c), (25) 
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where the vector V is external and, in the MSSM, stands for one of the neutral gauge bosons 
Z, 7. The coupling y t is the relevant Yukawa coupling proportional to m t and reads at lowest 
order 

= g m t 

Vt v^Mw/sin/?' 1 ' 

where g is the SU(2) gauge coupling. The various chiral fields are as follows. The SU(2) 

singlet u contains as physical fields the right component of top and its SUSY partner £r. 

The SU(2) doublet Q contains two chiral fields u, d, whose components are the left handed 

top and bottom, with their superpartners. The SU(2) doublet H u is one of the two Higgs 

doublets of the MSSM and contains two chiral fields H + , H° 











Q= \ 




H u =\ 






id) 







(27) 



Since V is external, the calculation can be entirely done in a gaugeless limit of the MSSM, 
precisely as in [22] where the subleading correction 0{G^mf) to the Z — > bb vertex and p 
parameter is computed in a gaugeless limit of the Standard Model. 

The Sudakov correction comes from the diagram proportional to g y\ shown in Fig. (j2J 
This is similar to the previous T v ^, but differs in the chiral structure because of the 




3>(P) 



FIG. 2: Diagram contributing the y\ Yukawa terms 

exchange of a $ field, instead of a vector V. The D-algebra details can be worked out 
and the expressions we derived for the universal correction are still valid apart from trivial 
changes in the couplings and, most important, the swap Cn <-> C21. This will be the 
technical reason for the correction being linear and not quadratic in the logarithm of the 
center of mass energy. 
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To perform the actual calculation, we should evaluate a few diagrams according to specific 
final state and exchanged virtual $ fields. Indeed, due to the explicit SU(2) structure of 
the Higgs field it is not possible to join together the photon and the Z. However, instead 
of considering the separate corrections to the 7 and Z vertices it is more instructive to 
relate the single logarithm appearing in c Y to the external $ self energy divergences. This 
is possible due to the following remarks. 

1. The linear logarithm and the divergence A of the vertex proportional to y\ appear in 
the fixed combination A — log(s/M 2 ). This follows from the actual expression of the 
diagram in Fig. (j2J). It is enough to substitute C\\ — ► C12 in Eq. (|2*2*j) and to combine 
the logarithm coming from the asymptotic expansion of C\i with the divergence of the 
B Q function from the first term of Eq. (|12|). 

2. The divergence in the vertex is canceled by the external $ self energies, also propor- 
tional to y^. This follows from the fact that the divergence of the sum of diagrams 
shown in Fig. Q and dictated by the SU(2) invariant e tensor is proportional to the 
sum of charges of the three chiral fields $1,2,3- Thus, cancellation results from the 
Yukawa interaction being a neutral SU(2) singlet, as it must. 




FIG. 3: Cancellation of the divergence in the terms proportional to y|. It is given by the sum of 
the above Yukawa diagrams where $1,2,3 can be H u , Q, or u as allowed by the SU(2) invariant 
coupling u Q ■ H u . 
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The above observations link the calculation of the logarithmic Yukawa correction to the 
external propagators which are are quite simple and independent from the off-shell vector V. 
Indeed, the coefficient of the logarithm in the vertex must be equal to the coefficient of the 
divergence in the self-energy. The correction to the kinetic term J d A 9 $$ in the effective 
action is well known and reads 

7l> 



where 



$) = y 2 I ^ A(p) d*9 S(-p, 9) S(p, 9), (28) 



A (P) = TF^ A + finite - ( 29 ) 



In addition, counting the possible diagrams and the color factors, we have the following 
multiplicity factors for the superfield propagators 

F(u,u*) : 2, (30) 

r(<5 Q ,Q ( g) : 8 a p, (31) 

T(H ua ,H uP ) : N c 5 aPl (N c = 3). (32) 

In conclusion, the relative Yukawa correction to the various processes one can consider is 
obtained in terms of the common coefficient 

q 2 m? , s a fh 2 , s 

FT lo S —2 = -^TTJT log =5> ( 33 ) 



2-16n 2 M 2 v b m t 2 8ns 2 w M 2 , 



m 



as in the following list 

c Y (u R ,u R ) = c Y (u R ,u* R ) = 2c, (34) 

c Y (u L ,u L ) = c Y (u L ,u* L ) = c, (35) 

c Y K,rf L ) = c Y (d L ,d* L ) = c, (36) 

c Y (#+ F") = c Y (#+, #-) = 3c. (37) 

If we now add the y 2 term — y b d Q- H d (where the second Higgs doublet is H d = (H%, H d ~) T ) 
we recover the full Yukawa correction Eq. (jSJ) for final fermions or sfermions. In particular, 
we proved SSR for this kind of correction with almost no calculations. 

For final charged Higgs bosons, we must recall that the physical H + field is defined as the 
combination H + = cos(3H^ + sin/3if7 . Therefore, for H + H~ production, we have to scale 
m t — > m t cos 2 (3 and m h — > nib sin 2 (3 in the above expressions obtaining precisely Eq. Q. 
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V. ANGULAR DEPENDENT SUDAKOV CORRECTION 



We now turn to the final piece of the calculation, the angular corrections coming 
from gauge boson exchange in boxes. The effective four point vertex is proportional to 
J d 4 6 $ 1 $2 < ^ ) 3 C &4) where $j are the chiral and antichiral fields required to describe the initial 
and final states of the scattering process. For instance, in four fermion scattering at tree 
level, the sum of 7 and Z exchange leads to the Dirac structure 7 M £3> 7 M . For right handed 
fermions in two spinor notation, this is proportional to <7 M ® <r M that can written in the form 
(1P11P4) (1P21P3) by a Fierz transformation. This is precisely the fully fermionic component of 
the above superspace integral. 

The angular dependent Sudakov correction comes from the two diagrams drawn in Fig. 
with their detailed chiral structure and from the two diagrams involving the l/ 2 $$ vertex 
shown in Fig. (J3J). For simplicity, we are assuming that the initial and final states do not 
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FIG. 4: Direct and crossed box diagrams giving the angular Sudakov correction. 




FIG. 5: Box diagrams with one four leg vertex V 2 &&. 
belong to the same (lepton or quark) family. 

nn 

For details about D-algebra manipulations, we refer to |lfil . I17l |. We just illustrate the 
general algorithm at one loop since box calculations with superfields are not usual. The 
various D 2 and D operators are moved away from their position by integration by parts in 
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8 space. They can be transferred to other internal lines or to external fields. When all the 



2 

D 2 and D on internal lines are on the same line, we obtain a zero result if there are less 



then four factors. If we have precisely D 2 D" , we can simply remove them. In particular, the 
D-algebra of Fig. (JHJ) is trivial. If we have more than four D, D factors, we can exploit the 
D-algebra {D a , D a } = ^i^aa to reduce their number. Useful relations are D^D^D^ = 0, 
D 2 D D 2 = 16D-D 2 and so on. After doing the D-algebra, we are left with a single d 8 
integration and a standard loop integral. 

The direct box diagram is reduced to Fig. (JBJ). As in the calculation of the universal 



K/v\/\/\yN 

D D.^ 4 D 

a B a P 

FIG. 6: Direct Box diagram: surviving term after D-algebra manipulations and neglecting mass 
suppressed terms. 

correction, we neglect systematically all mass suppressed contributions D 2 §, D $. Also, by 
definition _D$ = D$ = 0. Therefore, we can simplify D a D a §(p) ~ 2p aa §(p). The diagram 
is thus proportional to 2p^ ■ p^ Dq where Do is the four point basic PV function defined by 

I ^ D^Unr D(h m 

with 

D 1 = k 2 -M 2 ,D 2 = (k+p 1 ) 2 -M 2 , (39) 

D 3 = (k + Pl + p 2 f - M 2 , D i = (k + p 1 +p 2 + p 3 ) 2 - M 2 . (40) 

For simplicity and somewhat heuristically, we put the same mass M in all propagators. 
The final result will be independent from M at next-to-leading logarithmic level. We could 
of course, keep separate masses inside the loop and find the same result. We define the 
Mandelstam variables as s = (pi + p^) 2 , t = (p\ +p 2 ) 2 , u = {p\ +p 3 ) 2 . In the large s limit, 
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with fixed ratios t/s, u/s, the asymptotic expansion of the Do function is 

^-t^m^m)- (41) 

The relative value of the one loop correction with respect to the tree scattering amplitude 
is then 

n 2 / - » - t \ 

(42) 



direct 



^ 1 + log 2 



32tt 



M 2 



Ml 



To evaluate the crossed diagram it is convenient to perform D-algebra making all the 
D, D operators act on the fields $i and $4. This leads to the combination shown in 
Fig. (JZJ). Apart from trivial normalizations and couplings that we shall insert in the end, 



D <£, 



"P + 1/2 



(a) 



(b) 



FIG. 7: Crossed Box diagram: surviving contributions after D-algebra manipulations and neglect- 
ing mass suppressed terms. 



the diagram (a) gives a Co function |2l( due to the canceled propagator between fields 2 and 
3. The same is true for the two diagrams in Fig. (jSJ). The Co function and its expansion for 
s — (Pi + P2) ^ Pi, m i are 

C (pi,P2) = [ ^JTo 2\tn I \i 1 2\rn I 1 \ 2T ~ 7T lo S 2 — f ( 43 ) 



m 2 (A; 2 — m\)({k + pi) 2 — m|)((A; + pi + P2) — ml) 2s 



m 



Keeping only the terms growing with energy we obtain 



(a) + Fig. © = llog 2 ^. (44) 



Diagram (b) requires the expansion of the Du functions defined by 

d 4 k k» 



m 2 D 1 D 2 D^D i 



i=l 

These are given asymptotically by 



B "~iH 2 ^- 2iog2 ^)- < 46 > 

d "~-^Tp- (48) 
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As a final step, we must evaluate 

(b) = ~ J d'6 D a ^ $ 2 $ 3 D^ A {a^f a {D xl Pl + D 12 p 3 + D 13 p 3 ) M , (49) 

where Du = Du + D . We already remarked that (ft) aa D a Q(p) and its conjugate are mass 
suppressed. Thus, we manipulate each of the terms in Eq. (j4Uj) in order to have a D (or 
D) operator on each derived field (i.e. fields accompanied by a factor of the associated 
momentum). 

The term ~ p\ is already in this form and does not contribute. The ~ p 3 term is treated 
as follows 

" d A e D a ^ $ 2 $ 3 D^S 3 f a = (50) 

: - f d*6 D & D a ^ $ 2 $ 3 ^SX* + • • • = 4P1 • P3 [ d A 6 $! $ 2 $ 3 $4, 



where we used integration by parts of the D, D operators, because of momentum conserva- 
tion. In a similar way we obtain for the ~ p 2 term 



/ 



d 4 9 D a ^ $ 2 $ 3 D^ 4 (i> 2 f a = -Ap 2 ■ p A J d 4 9 $i $ 2 $ 3 $4- (51) 



Hence, 



W = l2 B (5 12 - Cl3 ) = l(log^ + 21og^). (52) 

(a) + Fig. © + (6) = 1 (log 2 ^ + log 2 ^) . (53) 

Reinserting factors of 2 and couplings, we find the full crossed contribution relative to the 
Born amplitude 

CL^K^ + log^). (54) 
Summing it with the direct term we obtain 

r ang _ ang ang _ 9p ( , 2 J_ _ , 2 NLO 9p , _j>_ , * 

C " direct + C crossed " 327r2 ^°g M 2 M 2 / ~~ 167T 2 Mz U ^ 

where we neglect subleading terms ~ log 2 (— t/s) or ~ log 2 (— u/s) not growing with energy. 
Also, in the single logarithm, we have chosen arbitrarily the scale Mz which is not fixed at 
this order in the logarithmic expansion. 

After the already discussed identification g — 2 (g'Y/2), we find perfect agreement with 
the calculation in components. Expanding J d A Q $ 1 $ 2 < ^ ) 3 <: &4 i n component fields, we recover 
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and 



SSR for the angular Sudakov correction. In addition, we have proved that SSR can be 
extended to the amplitudes for other four fermion or sfermion scattering processes, like for 
instance the fully scalar scattering f f* — > ff*. This is a pleasant result that can be checked 
by explicit component analysis [23]. 

As a technical remark, we notice that from the above expressions, it is straightforward 
to check that there are no box contributions of Yukawa type giving asymptotic Sudakov 
logarithms. The diagrams proportional to yf with exchange of virtual $ superparticles can 
be calculated with the same techniques and all the terms growing with the energy cancel. 



VI. CONCLUSIONS 



We have shown that superfield perturbation theory can be applied to the calculation of 
Sudakov electroweak corrections in the MSSM. The precise physical origin of these correc- 
tions permits to identify the supergraphs that produce logarithmically enhanced effects. For 
simplicity, we have analyzed the neutral, essentially abelian, contributions. The calculation 
is made easy by the physical nature of Sudakov corrections which are well definite, finite and 
gauge invariant. With a single supergraph, we simultaneously calculate the correction for 
several possible external states. In this framework, it is therefore quite straightforward to 
prove the existence of supersymmetric relations among the Sudakov corrections to various 
processes. 

It is clear that many extensions are possible, in particular to processes relevant to LHC. 
For instance, single top quark production via the parton processes hu — > td, ud —>■ tb, 
bg — > tW~ , and bg — > tH~ receives large SUSY electroweak corrections as discussed in [^] 
at the Sudakov level. The superfield approach can prove in such cases supersymmetric 
Sudakov relations with processes of production of single stop squark that have similar large 
corrections. 

As a final comment, we emphasize the remarkable simplicity of the Yukawa correction in 
superfield language and the possibility of computing it in a gaugeless limit of the MSSM. 
Indeed, it should be possible to evaluate the phenomenologically important two loop contri- 
bution in a rather compact way, at least in the Sudakov approximation, by combining recent 



results on multi-loop asymptotic expansions 
is in progress on such topic. 



24j with superfield perturbation theory. Work 
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